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Abstract — In  this  paper,  we  study  the  plane  wave  scattering 
from  perfectly  electric  conducting  (PEC)  bodies  of  revolution 
(BOR)  with  tip  singularities.  It  is  known  that  solutions  to 
surface  integral  equations  such  as  magnetic,  electric,  and 
combined  field  integral  equations  (MFIE,  EFEE,  and  CFIE, 
respectively)  are  singular  near  the  tips.  Consequently,  the 
convergence  of  method  of  moments  (MoM)  based  on  those 
surface  integral  equations  is  not  optimal  or  guaranteed.  By 
using  appropriate  graded  meshes,  one  can  retain  the  optimal 
convergence  rate  in  MoM. 

Keywords — Tip  Singularity,  Graded  Mesh,  Optimal 
Convergence,  Integral  Equations,  Method  of  Moments,  Radar 
Cross  Section. 


I.  INTRODUCTION 

Computational  electromagnetics  (CEM)  technology  has 
made  tremendous  progress  in  the  last  decade  due 
largely  to  the  advancement  of  fast  solver  and  high 
performance  computing  technology.  Consequently,  CEM 
tools  are  being  applied  to  ever  more  complex  problems. 
Even  though  CEM  tools  still  rely  on  radar  cross  section 
(RCS)  measurements  for  validation,  the  measurement 
community  is  increasingly  relying  on  CEM  tools,  especially 
those  based  on  the  method  of  moments  (MoM),  to  validate 
their  measured  data  to  minimize  measurement  uncertainty. 
To  provide  prediction  data  for  measurement  validation,  one 
typically  needs  to  compute  the  RCS  from  2 to  18  GHz  using 
very  fine  grids  to  ensure  solution  convergence.  This 
presents  quite  a challenge  for  MoM  codes  even  for 
canonical  targets  of  moderate  sizes,  especially  if 
computation  is  required  at  every  10  MHz  and  every  0.1 
degree.  Instead  of  using  a uniform  mesh,  one  would  like  to 
use  a non-uniform  mesh  that  is  denser  near  the  singularities 
and  coarser  elsewhere  to  minimize  the  number  of 
unknowns.  In  this  paper,  we  investigate  the  choices  of  non- 
uniform  mesh  that  give  fast  converging  solutions  for  bodies 
of  revolution  (BOR).  The  theorem  that  defines  the 
constraint  for  the  graded  mesh  will  be  discussed  and 


followed  by  numerical  examples,  particularly  those  of  a 10- 
foot  ogive  and  a 10-inch  ogive  with  gap. 


II.  Solutions  near  geometry  singularities 


It  is  known  that  the  solution  of  the  scattering  problem  by  the 
perfectly  conducting  ogive  is  singular  due  to  the  ogive  tips. 
In  fact,  if  r is  the  distance  from  the  tip,  the  solution  near  the 
tip  behaves  as 

where  and  //,  depend  on  the  angle  of  the  tip  (see,  for 

example,  [1]  and  [2]).  If  piecewise  polynomial  basis 
functions  defined  on  a uniform  mesh  are  used  to 
approximate  the  solutions,  the  convergence  of  MoM  is  not 
optimal  due  to  the  singular  behavior  of  the  solutions  near 
the  ogive  tips.  To  retain  the  optimal  convergence  rate,  one 
can  either  include  the  singular  basis  functions  in  the 
approximation  or  discretize  the  ogive  with  a graded  mesh.  It 
is  easier  to  construct  graded  meshes  and  apply  them  to  the 
existing  MoM  codes.  Here  we  apply  graded  meshes  to 
Cicero  [1]  which  is  a MoM  computer  code  for  bodies  of 
resolution  (BOR)  (see  [4],  [5],  [6],  and  references  therein). 
We  will  need  the  following  result  in  approximating  the 

singular  function  r~“  on  [0,1] , 0 < CU<  1 (see  [7]): 
Theorem  1:  Let  p be  such  that  (X  < 1/p . Define  the 
following  partition  T*  of  [0, 1] 


r.={Hny,  I = 0,1, 2, 3,...,  n, 
where  q = and  is  called  the  grading  exponent  of  . 
Let  S (t’  ) be  the  set  of  functions  that  are  constant  on  each 
subinterval  Then 

inf  (|r'“-^(r)||  =0(l/n). 


In  fact,  if  ^(r)  e 5 (f®  ) is  such  that 


ar)  = rr 


in  [r. , z;.^,  ],  i = 0, 1, 2,. . .,  n - 1 , then 
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The  norm  of  / e L''[0,l]is  ||/|Lr[o,,]  = '^j[|/|  dx  . In 

engineering  application,  I?  -norm  (i.e„  p = 2)  is  usually 

used.  In  this  case,  Theorem  1 states  that,  when  r”“ , for 
0 < r < 1 , and  OC<  )4,  is  approximated  by  step  functions, 

the  best  approximation  in  l3  -norm  is  obtained  when  these 
step  functions  are  defined  on  a graded  grid  with  grading 
exponent  q = . Note  that  the  above  theorem  can  be 

applied  to  singularities  of  arbitrary  order  such  as  vertices, 
edges,  comers,  etc.  The  singularity  order  (X  of  field 
solutions  near  the  ogive  tips  can  be  approximated  by  that  of 
the  field  near  the  tip  of  the  cone  whose  interior  angle  is  the 
same  as  the  angle  of  the  ogive  tip. 


III.  Graded  mesh  and  Convergence 


Let  S be  the  surface  of  the  PEC  ogive  which  is 
parameterized  by  arc-length 

(/,  <t>)  ^ (/?(/),  (f>,  Z{1)),  I e [0,  L],  <!> e [0,  In) 

where  L is  the  total  arc-length  of  the  generating  curve  and 
the  axis  of  rotation  is  along  the  z.  direction.  The  tips  occur 
at  /=0  and  l — h with/?(0)  = /?(L)  = 0 . The 

components  J,  and  J ^ of  the  surface  current  J = tixH 

behave  like  1'"  and  — near  the  tips,  where 
0 < Cr  < 1 . A graded  mesh  is  constructed  as  follows.  First, 
we  divide  [0,  L]  into  three  subintervals  [0,  £] , [f , L — £] , 
and  {L  — E,L\  where  £<LI4  and  is  a rational.  Let 
fi(x)  be  a twice  differentiable  function  defined  as 


£{x/£y 


b{x) 


for  xG[0,f], 
for  xe  {e,L-e], 

for  xe  [L-£,L], 


where  q = , 0 < p <\IOC , and  b{x)  is  a “connecting 

function”  which  is  monotonically  increasing  and  has  two 
continuous  derivatives.  (Such  a function  is  constructed  from 
a perfect  spline  in  [8]).  Then  the  nodes  in  the  graded  mesh 

is  defined  as 

l^=fi{x.),  x.  = L{i/N),  i = 


where  X.  ’s  represent  a uniform  mesh  and  /;  ’s  are  the 
mapped  points  in  the  graded  mesh.  Note  that  £ needs  to  be 
chosen  so  that  £ and  h — £ coincide  with  one  of  the  X-  ’s. 
It  is  common  in  MoM  codes  to  approximate  the  fields  with 
pulse  functions  (or  piecewise  constant),  i.e.  S(T^)  . Let 

: L'’([0,  L])  — > 5(r^)  be  the  orthogonal  projection, 
that  is. 


= Vv„eS(T^). 

Then  Galerkin  approximation  problem  is  to  find 
S (T^  ) such  that 

PAu„  = PJ, 

n n n*'  ’ 

where  A is  the  integral  operator  defined  as  in  the  combined 
field  integral  equation  (CFIE)  and  f is  the  given  incident 
field.  We  assume  that  the  integral  equation  Au  = / has  a 
unique  solution.  It  can  be  shown  that  Galerkin 
approximation  scheme  is  stable  [9],  that  is, 

IIPAv  |L  >Clivll, 

II  « «II/,-(0,2t,/,'’|0,L|)  “ II  "llz.-(0,2;r.z.'’|0,/.|)  ’ 


for  all  V„  e S {th ) and  some  C > 0 independent  of  . 
Consequently,  we  obtain  the  error  estimate  [9] 


u-u. 


"llL^0,2/r,L'’|0.Z.l) 


< C inf  M - V 


''llz,-(0,2/r,L'’|0./.|)  ’ 


where  u is  the  solution  of  the  continuous  problem  and 
S(T^)  is  the  approximating  solution.  This  implies 
that  convergence  rate  for  the  graded  mesh  is  optimal  (for 
L’’  -norm).  In  other  words,  the  approximating  solutions 
converges  to  u at  the  rate 

elements  in  It  is  possible  to  show  that  the 

optimality  holds  for  a more  familiar  weighted  Sobolev 
space,  L^([0,  L]) , whose  norm  is  defined  as 


Hence,  if  7,  ^ and  ^ are  the  approximate  solutions  of 
7,  and  7^ , respectively,  then  the  following  error  estimates 
also  hold 


7,-7,„  , , <C/n, 

II  ' Z.-(0,2^,Z.- |0,L|) 


'^«>.''IL3(0,2;r.4|0.Ll) 
where  (X  is  the  order  of  singularity. 


< C/n, 


rV.  NUMERICAL  RESULTS 

In  this  section,  we  illustrate  the  benefits  of  using  properly 
graded  meshes  with  Cicero  code  in  computing  RCS  of  a 10- 
foot  ogive  and  a 10-inch  ogive  with  gap.  Cicero  is  a MoM 
code  with  pulse  basis  and  testing  functions.  In  each  of  the 
following  figures,  the  numbers  in  the  legend  are  the  number 
of  sampling  points  per  wavelength  (ppw)  used.  For  graded 
meshes,  the  “ppw”  means  that  the  total  number  of  grid 
points  are  the  same  as  that  of  the  uniform  mesh  with  the 
“ppw”.  All  results  are  computed  using  CFIE. 

A.  The  10-foot  ogive 

The  10-foot  ogive  is  10-foot  long  from  tip  to  tip  and  1-foot 
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wide  at  the  waist.  The  tips  of  the  ogive  are  at  z = ±1 .524 
meters.  The  meshes  are  either  uniform  (q  = l)  or  graded 
with  the  grading  exponents  q = 2.  We  first  compute  the 
electric  current  components  and  7^  at  0.5  GHz.  In 
Figure  land  Figure  2,  we  plot  | /,  | , | 7^  | for  ^-polarized 
incident  field  at  0'  = 0°  and  6'  = 20° , respectively. 


d^Iarizdd  incident  field  at  6 


\ 

\ 

\ 

/ 

\ 
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Figure  1:  Magnitude  7^  and  7^  for  ^ -polarized  incidence  field  at 
^ = 0“  . Tips  of  the  10-foot  ogive  are  at  z=l  .524  m and  z=-l  .524  m. 


Due  to  numerical  limitations  in  Cicero  such  as  piecewise 
approximations  to  the  roof-top  basis  functions,  we  can  only 
observe  that  7,  and  7^  and  their  derivatives  tend  to 

infinity  at  the  tips  instead  of  become  infinity  as  expected  in 
[2].  In  any  case,  these  singularities  cause  slow  convergence 
in  MoM  using  uniform  meshes. 


z (meters) 


Figure  2:  Magnitude  7,  and  7^  for  0 -polarized  incidence  field  at 
^=20". 


We  assume  that  the  electric  currents  computed  with  640 
points  per  wavelength  are  “exact”  and  plot  the  relative 
errors  in  Figure  3 and  Figure  4 


fK,m 


X''^(0,L) 


Jy' 

L"-(0,L) 

, n = l,2,...,5 


where  V is  either  t or  (f>,  is  the  mesh  size 
corresponding  to  the  number  of  points  per  wavelength  pww 
= 2"xl0,  and  is  the  mesh  size  for  n = 6 (or  640 

points  per  wavelength).  It  is  observed  that  the  errors  in 
graded  meshes  are  smaller  than  their  uniform  counterparts 
for  sufficiently  large  ppw  and  decrease  at  a faster  rate. 


% 
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graded  meshes  are  smaller  and  decrease  faster  than  those  in  uniform 
meshes. 


6-pol  incic5ent  field  with  angle  = 20 


Points  perwa\«tenglh  (2"  x 10) 


Points  per  wavelength  {2"  x 1 0) 

Figure  4:  Logarithmic  relative  error  for  0 — 20*^  in 

graded  meshes  arc  smaller  and  decrease  faster  than  those  in  uniform 
meshes. 

Next  we  compare  the  RCS  results  for  different  meshes.  It 
can  be  seen  from  the  below  figures  that  the  graded  mesh 
performs  better  than  the  uniform  one  with  the  same  number 
of  unknowns,  that  is,  RCS  results  of  the  graded  meshes 
converge  at  a faster  rate.  In  Figure  5.  RCS  at  0.5  GHz  for 
69  polarization  computed  with  the  640  ppw-uniform  mesh 
is  viewed  as  the  “exact”  result.  The  zero-degree  angle 
corresponds  to  the  tip  scattering  direction  while  the  90- 
degree  angle  corresponds  to  the  broadside  scattering.  We 
observe  that  the  result  computed  with  the  graded  mesh  at  60 
ppw  already  overlaps  with  the  “exact”  solution  while  the 
RCS  curve  computed  with  the  uniform  mesh  at  60  ppw  has 

not  converged,  especially  near  9 = 20“ . In  fact,  it  requires 
at  least  160  ppw  in  a uniform  mesh  to  yield  the  same 
accuracy  as  the  60  ppw  in  a graded  mesh. 


figure  5:  RCS  of  the  10-foot  ogive  at  0.5  GHz  (99  polarization).  The 
“exact”  solution  is  represented  by  the  640  ppw-uniform  mesh  (solid  line). 
At  60  ppw,  the  uniform-mesh  solution  has  not  converged  (dotted  line) 
while  the  graded-mesh  solution  (dash-dotted  line)  overlaps  with  the 
“exact"  solution. 
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At  higher  frequencies,  tip  singularity  becomes  more 
problematic.  In  Figure  6,  RCS  at  5 GHz  for  the  10-foot 
ogive  is  computed  near  the  nose-on  (grazing)  angular  region 
by  using  uniform  meshes  of  different  grid  densities.  The 
curves  oscillate  near  -69.5  dBsm. 


Figure  6:  Uniform  mesh  - RCS  of  the  10-foot  ogive  at  5 GHz  near  the 
grazing  angle  ( ^ polarization). 


However,  RCS  at  5 GHz  for  the  10-foot  ogive  computed 
with  graded  meshes  approaches  monotonically  to  a 
converged  solution  as  seen  in  Figure  7. 


Figure  7:  Graded  mesh  - RCS  of  the  10-foot  ogive  at  5 GHz  near  the 
grazing  angle  ( ({xj)  polarization). 


angle  of  22.62  degree,  and  with  a maximum  radius  of  1 inch 
at  the  middle.  The  generating  curve  for  this  ogive  is  part  of 
a circular  arc  with  a 13  inch  radius.  A small  rectangular 
groove  cut  out  around  the  middle  of  the  ogive.  The 
circumferential  groove  is  0.25  inches  wide  by  0.25  inches 
deep.  The  bottom  of  the  groove  forms  a ring  0.75  inches  in 
radius. 


Figure  8:  Ogive  with  gap. 


Here,  the  solutions  to  the  integral  equations  have  both  edge 
and  tip  singularities  due  to  the  groove  and  the  ogive  tips.  As 
in  the  case  of  the  10-foot  ogive,  we  design  appropriate 
graded  meshes  to  improve  the  convergence  rate  in  MoM 
solutions.  There  are  tlu-ee  comers  in  the  generating  arc  and 
each  has  different  angles.  Thus,  we  construct  a graded  mesh 

with  three  grading  exponents  (^p92>9.s)*  example, 
the  distribution  of  points  in  a graded  mesh  with 
( ) “ ( ^’  2, 2)  is  plotted  in  Figure  9. 

0.025 
002 

_ 0.015 

E, 
a 

0.01 
0.005 
0 

Figure  9:  Point  distribution  of  a graded  mesh  on  the  generating  arc  of 
ogive  with  gap. 


Graded  mssh  for  the  ogive  with  gap  ([q,  .Pj  ,q  jj  = [2 ,2 ,2]) 


B.  The  10-inch  ogive  with  gap 

To  examine  further  the  advantage  of  graded  meshes,  we 
also  consider  the  ogive  with  gap  (see  Figure  8),  one  of  the 
test  targets  proposed  by  Electromagnetic  Code  Consortium 
[10],  [11].  The  ogive  is  10  inches  long,  subtending  a half 


In  Figure  10  and  Figure  11,  RCS  curves  at  2 GHz  are 
plotted.  Uniform-mesh  solutions  converge  much  slower 
than  graded-mesh  solutions,  especially  near  the  grazing 
angular  region.  We  see  that  at  80  ppw,  RCS  of  the  graded 
mesh  already  converges  while  that  of  uniform  mesh  does 
not. 


98 


ACES  JOURNAL,  VOL.  19,  NO.  2,  JULY  2004 


Unnorni  muti  - og>v«  with  gap  (2GHz).  PP 


Figure  10:  Uniform  mesh  - RCS  of  ogive  with  gap  at  2 GHz 
( polarization). 


I 


'o  ie903040SOe070M90 

Figure  11;  Graded  mesh  - RCS  of  ogive  with  gap  at  2 GHz 
( (fx/)  polarization). 

Furthermore,  the  RCS  differences  at  ^ = 0°  (the  grazing 
angle)  decrease  for  graded  meshes  while  oscillate  for  the 
uniform  meshes  as  seen  in  Figure  12  and  Figure  13. 

V.  CONCLUSIONS 

In  this  paper,  we  present  a construction  of  graded  meshes 
that  enable  a faster  convergence  using  MoM  for  BOR 
targets  with  tip  singularities.  Numerical  results  are  given  for 
the  PEC  ogive  with  and  without  gap  using  CFIE.  These 
preliminary  results  show  that  faster  convergence  can  be 
achieved  if  one  chooses  a graded  mesh  using  the  technique 
outlined  in  this  paper.  We  have  also  observed  the  similar 
improvement  in  convergence  for  other  types  of  integral 
equations.  The  technique  can  be  easily  generalized  to  non- 
BOR  targets  which  have  tip  and  edge  singularities.  This  will 
be  reported  in  the  future. 


Figure  12:  Uniform  mesh  at  2 GHz  - RCS  errors  in  ogive  with  gap  at 
^ = 0", 


Figure  13:  Graded  mesh  at  2 GHz  - RCS  errors  in  ogive  with  gap  at 

e=Q\ 
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